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ABSTRACT 


This program has been concerned with the study and improvement of 
the short-term spectral purity and the long-term frequency stability 
of stable laser oscillators. 

We have successfully measured the short term frequency phase 
fluctuations in high quality He-We gas lasers^ and have measured^ for 
the first time, the ultimate limit on laser spectraJ. width caused by 
quantum phase noise, as predicted by the Sc hawlow- Townes formula. 
Verification of this formula, one of the most basic in laser theory, 
has only recently come within the range of experimental feasibility. 

■ An additional task of this program was to study and, if indicated, 
to demonstrate experimentally a new absolute laser frequency stabili- 
zation method. The study has been completed and experimental verifi- 
cation is underway, with many problems solved. 

Finally we have suggested an important new mechanism to account 
for the frequency fluctuations in He -Re 6328 X lasers, and presented 
both experimental and theoretical support for this mechanism. 




I. -COMPLETED PEEQUENCY FLUCTUATION MEASUEEMENTS 


This section of the report summarizes laser frequency fluctuation 
and. quantum phase noise measurements that were completed and published 
during the period of this contract. 

Even in a very high quality, stable, single-frequency laser oscil- 
lator, the instantaneous laser phase or frequency will retain some noise 
fluctuations, due to two different sources. First of all, and usually 
most important in practical situations, are various external disturbances, 
by which we mean such things as microphonics, structural vibrations, 
plasma oscillations and, other disturbances in the laser. In addition, 
as a fundamental or ultimate limiting source of -noise, there are quantum 
phase or frequency fluctuations that arise from the basic process of 
atomic spontaneous emission in the laser medium. The magnitude of the 
oscillator spectral broadening due to these ultimate quantum noise sources 
is given by the so-called Schawlow- Townes formula, indicated at the bottom 
of Fig. 1. 
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FIG. 1--Frequeney fluctuation sources. 
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One of our primary purposes in this work has been to observe these 
quanta noise fluctuations in a stable laser oscillator^ and thereby to 
verify this Schawlow- Townes formula. We have achieved this objective; 
and along the way have also learned some important things about the 
apparent primary source of external disturbances in the 6328 % He-lfe 
laser in particular. 

Figure 2 is a schematic showing the essential features of the 
experimental apparatus for these measurements.’ -As indicated, we hetero- 
dyned together two stable single-frequency high-quality He-Ne lasers, ■ 
which are basically Spectra-Physics Model 119 lasers, shorn of their 
regular cabinetry and equipped with special power supplies and control 
circuits. As indicated in the upper sketch, the frequency of laser 
number 1 is set very close to the edge of the laser oscillation range, 
in order to obtain stable operation of this laser at very low power 
level, thereby enhancing the size of the quantum noise disturbances 
because of the l/p dependence in the Schawlow-Townes formula. 

Two sets of experiments were carried out. The first of these used 
a beat frequency of 30 MHz and commercially available electronics . This 
work was reported in Phys. Rev. -Letters 1? (April 19^8), "Observation 
of Quantum Phase Noise in a Laser Oscillator"; Appendix I of this report. 
The second set enployed a U.5 MHz beat frequency and was reported in 
J. Appl. Phys. 40, (February 1969), "Further Measurements of Quantum 
Phase Noise in a He-Ne Laser"; Appendix II of this report. We will sum- 
marize the results of these measurements very briefly here. 

The two lasers are heterod3nied together on a photomultiplier -tube, 
with ajnominal four-and-a-half MHz or 30 MHz difference frequency, and 
the resulting beat note is fed into a 4.5 MHz or 30 MHz RF discriminator. 

- 2 - 



EXPERIMENTAL APPARATUS 



FIG. 2-“Schematic showing the essential features of experimental apparatus. 






Part of the output of this discriminator is fed bach through a slow auto- 
matic frequency control loop with a long time constant in order to keep 
the frequency of laser spaced away from laser by 4-l/2 MHz or 

30 MHz on the average. 

Due to the instantaneous frequency or phase fluctuations in each 
laser, there is a finite amount of frequency jitter in the beat frequency, 
with a spectral width and a peak deviation in the 10 KHz range. This 
frequency jitter is reflected in the audio output of the EP discriminator, 
which is taken to an audio wave analyzer, whose output vs .‘audio frequency 
is then recorded on an xy recorder. 

Figure 3 illustrates the general character of the discriminator audio 
output spectrum as recorded on the xy recorder. The effect of external 
disturbances is essentially to modulate the instantaneous frequency of 
the lasers, or of the beat note, in a G-aussian noise-like fashion; and 
the effect of this is to give a l/f audio output spectrum in the dis- 
criminator output, extending from the cutoff frequency of the AFC loop, 
v?hich is a small fraction of 1 Hz, upwards. By contrast, the effect of 
the quantum noise may be described as a random walk of the instantaneous 
phase of the laser oscillator. This has a different spectral character- 
istic, and appears as a white noise output in the audio output spectrum 
of the discriminator. Because of the l/p dependence of the quantum 
noise term, this quantum noise level increases with decreasing 'laser 
power, as indicated on the figure; while the external disturbances repre- 
sent frequency modulation of the laser by external sources and is inde- 
pendent of the laser's oscillation level. 
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FIGURE 3 



Figure k shows a typical experimental result^ indicating how well the 
experimental measurements agree with the theoretical expectations . Also 
in agreement with theory^ we observe that the external contribution portion 
of this’ curve does not change as the oscillation power level of the lasers 
is- changed; whereas the quantum contribution does indeed increase as the 
inverse of the laser power level,, We should emphasize here that in all of 
our experiments, the frequency fluctuations of the laser and the resulting 
spectral broadening of the laser output were dominated by the external 
contribution, which gave a total spectral spread for the lasers, and con- 
sequently for the beat note aJ-so, of about 10 KHz,. The quantum contribu- 
tion, on the other .hand, represented a net broadening contribution of only 
a few 100 of cycles or less . Nonetheless, this quantum contribution could 
be separated out and measured precisely because of its different spectral- 
appearance in the outputs of the discriminator , 

It should be noted that radio frequency discriminators are in general 
rather difficult devices, with far from ideal noise properties; and spurious 
results that appear very much like the expected quantum contribution results 
can be easily obtained simply as a result of the inherent noise properties 
of the discriminator. -Appropriate steps were taken to correct for these 
difficulties . 

By making a large series of measurements such as the typical example 
shown in Fig. k, with the laser operating at different power levels 

in each case, we were able to obbain the data on quantum linewidth contri- 
bution in Hz vs laser oscillation power in watts shown in Fig. 5* These 
represent the principle quantum noise experimental results of this report, 
and show quite clearly the expected l/p dependence predicted by the 
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FIGURE k 





FIGr. 5““Da'ta on quantim linewidth contribution in Hz vs laser oscillation 
power in watts . 




Schawlow-Townes formula over a range of nearly two orders of magnitude. 

The lower points represent the results of the first series of experiments 
using an Ef discriminator centered at 30 MHz. The upper curve represents 
our more recent results using the 4-1/2 MHz discriminatorj, together with 
a few experimental points obtained by measuring the beat note phase 
jitter using a sampling oscilloscope technique. Various laser tubes and 
mirrors were changed between the two series of measurements, and the 
approximate; factor of. '2. difference may be a good measure of the repro- 
ducibility of our results . 

The fuzzy line represents the prediction of the Schawlow-Townes 
theory, taking -into account the uncertainties in -how well we know the 
cold cavity bandwidth and other parameters of our laser. The excess 
noise ratio for our lasers was determined by carrying out a separate 
series of measurements of the excess amplitude fluctuations of the lasers 
at low powers, following the method of Fried and House. From such measure- 
ments the excess noise factor can be determined directly, and was found to 
have a value of approximately 3 » 

Note that the apparent noise contribution begins to increase sharply 
at about the signal power level at which the signal-to-noise ratio — that 
is the ratio of the beat signal to the photomultiplier shot noise — 
decreases below about 20 dB. We believe, however, that this represents 
excess noise behavior in the discriminator itself as the signal level 
drops below the necessary threshold level for good fm detection; rather 
than any real increase in the quantum noise fluctuations. 
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II > ON-GOING FEE(^UENCY H.UCTUATION ME3ASDEEMEHTS 


This section of the report sumMarizes some additional frequency 
fluctuation -measurements which were started and carried to a fairly 
advanced level, hut not fully completed during the period of this 
contract <, 

The- Schawlow“Townes formula .suggests that Af^^ may he increased 
to a detectable magnitude hy decreasing P and increasing Af 

' C3.V 

Decreasing P in de-excited lasers .generally leads to an increase in 

plasma noise which swamps out Af^^ ; therefore, rf excitation on rf 

quieting of the plasma is indicated. Increasing Af^^^ is accoH 5 )lished 

hy lowering the cavity Q with low reflection mirrors. Naturally the 

laser cannot he made to lase at all unless the threshold requirement is 

met; i,e», gain must he high enough to overcome mirror transmission, 

A laser line with one of the highest known gains is X = 3»508 |i in 

He-Xe, He-Xe jj, lasers produced here recently have shown gains 

in excess of 100 dB/meter, Further, we have found the He-Xe plasma 
I 

to he easily controllable over a wide range of Xe pressures. 

However, because of the very rapid clean-up of the Xe gas in the 
discharge, lasers of this type previously have been limited to an 
operating life measured in "hours. We have found. that an apparently 
completely success.ful solution- to this problem .is to heavily overfill 
the laser with Xe, and- then condense out the excess Xe in a tside arm 
cooled to near liquid nitrogen ten 5 )eraturd . The Xe actually solidifies 
in the side arm; .however, the vapor presstire of Xe at just above liquid 
nitrogen ten 5 »erature (77°l) is 19 to 20 microns of mercuiy; i,e,, 
optimum for lasing. 
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Measurements of the ~ 100 MHz wide 3*51 M- transition of the He-Xe 
laser show no clearly identifiable Lamb dip. The absence of the expected 
dip may be due to a large homogeneous collision broadening linewidth, of 
the order of the inhomogeneous doppler width, introduced by the high 
pressure chosen for maximum lifetime together with the fact that the Xe 
we have used so far is not isotopicly pvire. 

One of our objectives in this work is a direct measurement of M 

q. 

with He-Xe lasers. To get some idea of the experimental desirability of 
He-Xe over the common 6328 X He-He laser, the Schawlow-Townes prediction 
is, plotted for -several lasers in Fig. 6. 

Two 30 cm external -mirror He-Xe laser on invar bases have been 
fabricated and have undergone preliminary experimentation (Fig. 7). 

Gain per unit length at various partial pressures of Xe, optiminn mirror 
reflectivities, Brewster window materials, detector systems, and align- 
ment difficulties have been some of the problems under study. We have 
now seen output powers of 2-1/2 milliwatts at 3.508 p using 14gF„ or 

i 

CaFg Brewster windows and output mirror substrates coated hO^ reflective. 

Beats have been obtained from these lasers and locked at 200 KHz. 

The beat was very noisy and showed clear 60 Hz PM. Work is now in progress 
to quiet the electronics and improve the pressure stabilization on these 
lasers . 

Two dc excited Spectra-Physics ll9 lasers, formerly used for quantum 
noise measurements at 6328 2 were modified for 3»508 He-Xe operation 
(without retiarn paths to avoid cataphoresis effects in the Xe discharge) . ' 
Although these had good temperature-compensated invar cavities, they were 
extremely noisy both in amplitude and phase with power outputs of only 
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FIG- J--30 cm external mirror cavity. 
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20 to 30 microwatts. The bandwidth of the beat Vas larger than the 
1 MHz bandwidth of our detector (about 7 MHz). Conversation with 
Dr. Gamo of the University of California at Irvine who had a similar 
laser^ corroborated this result. Improvements in plasma tube design and 
specifically the inclusion of a gas retiirn path might reduce this plasma 
noise and bring the beat width within the detector bandwidth. 

Finally, a quartz ingot internal mirror laser is nearing completion 
(Fig. 8). This laser could serve as the prototype for internal mirror 
Cer-Vit lasers providing maximum mechanical and thermal stability. 

Cer-Vit is a new low-expansion glass ceramic with -an expansion coefficient 
adjustable through zero at 25 °C. 

Liquid nitrogen cooled InAs photo conductive detectors, can achieve 

-9 

20 dB signal-to-noise ratios at 10 watt power levels. As can be seen 

from Fig. 1, such a low -power level should render Af easily measxirable. . 

y. 

We are confident that measurement of Af over a wide range of P will 
be achieved, thereby checking the fundamental quantum noise theory. 

III. STUDY OF CAUSES OF FREQUENCY FLUCTUATIONS 

In Applied Physics Letters, vol. 13, No. 6, pp. 197 -I 99 (September 
1968) "Current Pushing of the Oscillation Frequency of a 6328 S. He-Ne 
Laser" (Appendix III), we reported our results in explaining what seems 
to be the primary source of the external contributions which still limit 
the overall laser spectral width of He-Ne 6328 % lasers . It should be 
noted that in our work, and also in most other work with 6328 % lasers 
that we have seen reported, it appears that as one improves the laser 
design, and improves the shielding of the laser from acoustic and micro- 
phonic disturbances, the laser spectral width decreases to something on 
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FIG. 8 — Ingot laser. 
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the order of 10 KHz; and then sin^jly does not get any better as further 
improveinents are made. No one seems^ for exan^le, to have been able at 
6328 % to achieve the beat note purity of 20 or 30 cycles obtained by 
Javan, Jaseja, and Townes in their "Wine Cellar Experiments" which were 
done at 1.15 M- We now believe that this is associated with an unusually 
high polarizability of the He-Ne discharge at 6328 X, caused by the 
presence of Neon 1-s metastable atoms, and by the resulting mod\ilation 
of the laser cavity frequency through turbulent disturbances in the laser 

I 

plasma. • ’ 

Figure 9(a) shows a sinple experiment that was done to investigate 
the current pushing in our lasers, that is the small shift M in 
instantaneous laser oscillation frequency caused by a small change Ai 
in the discharge current through the laser tube. Very small shifts 
caused by small changes in discharge current could be easily measured 
using the heterodyne technique and the rf discriminator. 

Figure 9(b) shows the experiment^ results for the current pushing 
figure, that is, the derivative of laser frequency vs laser discharge 
current, as a function of where the laser’s oscillation frequency is 
tuned within the atomic linewidth. Q?his quantity is fairly large, and 
absolutely independent • of vfhere the laser is tuned within its oscillation 
range. This indicates that the current pushing effect is not caused by 
any pushing or pulling effects associated with the lasing atomic transi- 
tion itself; but must be caused by a background change in the dielectric 
constant of the plasma discharge, having nothing to do with the laser 
transition. Similar measurements were made at 1.15 H also, as indicated 
on the bottom of -the figure, and it was found that the current pushing 
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figure in this case was of opposite sign, very much smaller, and not 
necessarily flat across the atomic line. We will concentrate here on 
explaining the 6328 S results . 

Various mechanisms for explaining this result were considered. For 
example, the change in free electron density with increasing current in 
the laser discharge is many orders of magnitude too small to account for 
the observed current pushing. What we believe to be the correct explana- 
tion for this effect is illustrated in Fig. 10, showing the relevant 
energy levels of the He and Ne atoms, notice that there is a collection 
of lowest excited We levels, the 1-s metastable levels, which build up a 
sizable population in the laser discharge. Furthermore, there are two 
relatively strong absorptive transitions at .633^ i and 64-02 £ from the 
most heavily populated of these metastable levels up to certain of the 
2p levels which are not directly connected with the laser transition. 

The numbers in parentheses indicate the oscillator strengths of these 
two transitions. Increasing the current through the laser increases 
the populations of the He l-s metastable levels; and the reactive 
effects of these excited atoms then change the effective dielectric 
constant of the cavity at the nearby laser transition. 

Figure 11 illustrates this point. The laser transition at 6328 S 
is far out on the tails of the dispersive ciirves of these two contri- 
butions. Wonethless, evaluation of the expected metastable density, the 
oscillator strength of these transitions, and the resultant pulling effects 
maKe it clear that these two transitions are the culprits, note that the 
weaker but much closer 6334 X transition contributes about two-thirds of" 
the effect while the stronger but more distant 64-02 S transition contri- 
butes the remaining one-third. These are absorptive transitions; and 
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FIG. 9 — (a) Heterodyne experiment for measuring small shifts in 

laser oscillation frequency caused by current modulation 
’ or other perturbations. 

(b) ' Shift in laser oscillation frequency with discharge 

current, df/dl , measured at both 6328 2. and' 1.15 

in the same lasers, as a function of laser frequency 

\ 

. tuning -within the atomic linewidth, for three different 
values of dc discharge current 
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FIGo 10 — He-Ne energy diagram. 


- 20 




absorptive transition "pushes" a cavity frequency away from it; and hence 
the shift in the cavity frequency is to shorter wavelength, that is a 
blue shift with increasing discharge current . An extensive theoretical 
discussion of this effect is included in Appendix IV. 

It is, therefore, our conclusion that the residual spectral broadening 
in high quality 6328 % lasers probably results from random modulation of 
the cavity resonance frequency caused by random noise and turbulance .in 
the laser discharge, acting through this Ne metastable mechmism. Note 
that in all of our experiments, care was taJfcen to see that no coherent or 
single -frequency plasma oscillations were present or detectable by any 
means in the laser apparatus. Even so, .there is still a substantial 
amount of low frequency l/f type of noise in the plasma discharge. No 
simple mechanism for eliminating these effects is apparent at present, 
and this effect must be considered in any future applications of the 6328 % 
laser, in particular as an ultra-stable frequency standard or laser com- 
munications source. 

IV. -ABSOLUTE PEEQUENCY STABILIZATION OP He-Xe - 3.51 p 'LASERS 

As part of the effort under this grant, we have also considered ; the 
long-term absolute frequency stabilization of the He-Xe 3. *51 P laser. We 
have proposed . that a frequency discriminant be obtained by frequency- 
modulating a portion of a primary laser oscillator's output and passing 
it through an lonsaturated and non-regenerative reference laser amplifier. 
The PM-to-AM conversion at the modulation frequency in the reference laser 
amplifier becomes zero only when the primary laser oscillator frequency 
exactly coincides with the atomic line center of the .reference amplifier. 
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Thus, the FM-to~AM conversion provides a frequency discriminant for a 
stabilization feedback loop. The basic stabilization scheme is shown 
in Fig. 12. 

This absolute stabilization scheme has the following desirable 
attributes : 

(a) The main laser output is completely free from any applied 

modulation sidebands or frequency since only the 

portion tapped off to the reference amplifier is modulated. 

(b) The reference laser amplifier operates in a linear, small- 
signal condition, thus avoiding any of the complexities which 
may be associated with saturation, "hole burning" effects, or 
other second-order effects » 

(c) The gas fill, pressure and pumping conditions of the reference 
laser amplifier can be adjusted for optimum long-term frequency 
stability and reproducibility; while these parameters in the 
primary laser oscillator can be adjusted for maxim'um output 
power, efficiency, or other considerations independent of 
frequency stability. 

This scheme should be particularly effective with the 3 »51 He-Xe 
laser because: 

(a) The He-Xe laser at 3*51 M- has high gain, which permits good 
sensitivity for the present scheme . 

(b) Xenon is a heavy molecule, and thus has a narrow doppler 
linewidtho 

(c) The He-Xe plasma is a cool plasma, and this also leads to a 
narrow dopper line. 
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(d) The He-Xe laser amplifier can be operated at very low pressure^ 
thus having minimal pressure shifts. 


(e) The liquid Ng pressure stabilization method used in our laser 
tubes should provide long-term stable operation of the lasers, 
with essentially constant Xe pressure. 

Three pure xenon amplifier tubes have already been constructed, one 
50 cm long and two 25 cm long. The inner bore diameters are 1 5 mm. 
Choosing this large bore diameter permits easy alignment of the laser 
beam through the amplifier tubes, while also lowering the gain factor to 
avoid saturation and regeneration effects. Each tube has a side arm with 
a reservoir bulb which is cooled down to liquid Ng temperature in opera- 
tion, so that the tube is operated at a pressure corresponding to the Xe 
vapor pressure at that temperature. All of the tubes are designed for 
dc operation, so that pumping conditions can be controlled. However, the 
two 25 cm tubes can also be operated under rf excitation. These tubes 
appear to be satisfactory for our e^jperiments . 

To date all components of this system except for the necessary phase 
sensitive detectors, power amplifiers, and associated feedback equipment 
have been designed and are under construction or on hand. 

Final reference amplifier gain measurements will begin as soon as the 
laser amplifier and oscillator can be set up in the vaixlt. The frequency 
modulation system will be tested when the LiNbO^ modulator is constructed. 
Finally, the whole stabilization scheme will be tested* and the short-term 

and long-term stability measured. From sensitivity calculations, a fre- 

3-IL~X2 

quency stabilization as good as one part in 10 is possible with this 

system. However, the difficulties introduced by various practical con- 
siderations, such as interference and imperfect components, can be deter- 
mined only with the acutal experimental system. 
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A theoretical discussion, of this system was presented in Applied 
Optics, vol. 5^ No. 10, October 1966 ^ and is included for reference in 
Appendix V. 

This completes the final report on Grant NGR 05(020) -23^. -All of 
the. aims of the original proposal have been met or exceeded and a firm 
foxmdation for ongoing and future research has been laid. 
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By heterodyning together two stable 6328-A He- 
Ne lasers, with one laser operating at very low 
power, quantum phase fluctuations caused by • 
spontaneous emission have been observed. Re- 
sults, although preliminary, seem in good agree- 
ment with the predictions of Schawlow-Townes 
and others.^ 

Neglecting amplitude fluctuations, the beat 
note between two laser oscillators may be writ- 
ten n= Vo cos (wof+^), where coq is the mean fre- 
quency and <p(t) the randomly varying phase. 

Since the iiBtantaneous beat frequency is Wo 

the power spectral density G^{/) of the 
quantity may be determined with an rf fre- 
quency discriminator centered at Wq, followed 
by an audio-wave analyzer (Fig. 1). With appro- 
priate instrumentation one can also measure the 
analytically related mean-square phase jitter 
(a^(t))= ([ 9 ?ft+T)-^(f)f) as a function of the 
time interval t. 

In real lasers the random phase variation (pit) 
includes an “external” contribution <pgit) due to 
acoustic noise, structural vibrations, and plas- 
ma disturbances, plus a usually much weaker 
contribiRion (pqii) due to quantum noise. In our 
ejcpertments the external disturbances occur pri- 
marily at low audio frequencies, with G^gif) 
~l/f‘ [see Fig. 2(a)] and (A<p/(t))~t^. The re- 
sulting beat-note power spectral density if) is 
Gaussian with a linewidth typically A/g *3.5 kHz 
in our apparatus, essentially independent of la- 



FIG. 1. Block diagram of system used to observe 
quantum phase fluctuaticms. 


ser power level. By contrast, the quantum con- 
tribution should have a white-noise spectrum 
G^qif) = 477 A/g and a phase jitter <Aq3g“(T)> 

= 27TAgT. The Schawlow-Townes prediction is 



nhfiAf ? 

cav 2 


where / is' the oscillation frequency,, A/gav 
“cold”-cavity bandwidth, P the laser-oscillation 
power level, and Nj and the i^per and lower 
level populations. Except at our lowest power 
levels, the quantum contribution to the total beat- 
note spectral density is considerably less than 
the Gaussian external contributions. The quan- 
tiun contributions are, however, separable by 



FIG. 2. (a) Typical discriminator noise output spec- 
trum as measured by audio-wave analyzer, showing 
1/f^ portion due to external disturbances and flat por- 
tion due to quantum noise (or, in some oases, to dis- 
criminator characteristics), (b) Quantum phase-noise 
linewidth contribution, as measured by flat discrimina- 
tor noise level, versus oscillation power level of laser 
i-j. Also shown are the Schawlow-Townes theoretical 
result assuming Wj/fNj— Wj) = 1 and the results of exper- 
imental phase-jitter measurements at a single fixed de- 
lay t= 167 nsec. 


ReproductkJii in whole or in part is permitted by the iwblisher 
for any purpose of the United States Government. 
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examining G^(/) at high enough frequencies, or 
at short enough times . 

The- experimental apparatus (Fig. 1) is basical- 
ly the same as reported earlier.* Laser has 
a relatively high-transmission output mirror 
(1.7%) to enlarge and thus enhance the 

quantum line broadening, A slow power -stabili- 
zation loop provides stable operation al low pow- 
er levels by piezotuning the cavity close to the 
edge of its oscillation range.* Laser Lj is 
locked 30 MHz from Lj by a slow automatic-fre- 
quency-control loop. The signal-to-noise ratio 
(S/N) of the beat-note depends upon signal power 
P from Lj and photocathode quantum efficiency 
in the usual way. Since these measurements nec- 
essarily involved low and decreasing S/ff, it was 
necessary to verify that measured noise line- 
width increases at low P did not simply rep- , 
resent equipment characteristics. This was 
checked following each measurement by operat- 
ing Li at a higher ou^ut level, where quantum 
fluctuations in Li should be negligible^' and in- 
serting a variable optical attenuator (~30 dB) in 
the output of Li to produce the same post-attenu- 
ator power output (and hence S/N conditions) in 
the succeeding apparatus. Comparison of the 
two measurements effectively determined the 
quantum contribution. 

The circled points in Fig. 2(b) show the mea- 
sured white-noise level of G^if), expressed as 
equivalent linewidth A/^,-. versus oscillation lev- 
el of Lj. The theoretical curve is the Schawlow- 
Townes formula* taking into account the uncer- 
tainty in cavity par^eters for L, but assuming 
Va/ (Va-gaV^/gi) = 1 . The offset between theory 
and experiment can be accounted for by assum- 
ing not unreasonable’ for this 

particular laser system. The square data points 
represent measurements of (Asp*(t)> at one fixed 
value T = 167 nsec, again converted to. equivalent 
quantum linewidth. As also in our earlier work, 
we are unable to resolve the factor of 2 differ- 
ence between discriminator and phase -jitter re- 
sults here, and must continue to attribute it to 
eiqperimental uncertainties or to some systemat- 
ic error in one of the measurement techniques. 
We discount the apparent rapid increase in aA 
below P = 2x 10“'* W since in this region the beat- 


note S/N decreases below 20 dB. At such low 
power levels, fluctuations in the power stabiliza- 
tion loop could result in an increase in observed 
noise because of nonlinear power dependence 
(l/P) of the quantum noise. 

Our useful measurement range is uncomfort- 
ably limited at present by the vanishing quantum 
contribution at higher values of P and by reduc- 
tion of the heterodyne S/N at lower values of P. 
However, we believe that the 1/P dependence ob- 
served between P=2x 10“* W andP=9xlO“*W 
most probably represents quantum phase fluctua- 
tions in laser L^. In future experiments more 
detailed study should be possible by using a high- 
.gain infrared laser transition, reducing the cavi- 
ty Q, increasing A/cav> thus greatly enhanc- 
ing the quantum noise contribution. 

Note added in proof . — Just as this report was 
completed we received the translation of a Rus- 
sian Letters journal reporting very similar ob- 
servations,^ although at substantially higher pow- 
er levels P and hence much lower values of A fq 
(-0. 1-1.0 Hz). Our only reservations concern- 
ing the Russian results have to do with the non- 
ideal discriminator characteristics mentioned 
earlier; i.e,, we find that the observed inherent 
noise output from our real (nonideal) discrimina- 
tor can have a flat spectrum that rises as ~1/P 
due simply to decreasing S/N in the rf bandwidth 
rather than to any real frequency fluctuations in 
the beat signal. The results of Ref. 4 imply the 
availability of a very nearly ideal rf discrimina- 
tor at = 2?r X 8.4 MHz. 


♦Work supported by Joint Services Electronics Pro- 
gram Contract No. Nonr-225(83), and by a National 
Science Foundation Traineeship for R. Arrathoon. 

*A. L. Sohawlow and C. H. Townes, Phys. Rev. 112. 
1940 (1958); for a more recent discussion, cf. M. Lax, 
in Physics of Quantum Electronics, edited by P. L. 
Kelley, B. Lax, and P. E. Tannenwald (McGraw-Hill 
Book Company, Inc., New York, 1966), p. 735. 

^A. E. Siegman, B. Daino, and K. R. Manes, IEEE J. 
Quantum Electron. QE-3 . 180 (1967). 

*F. T. Arecchi, A. Berne, A. Sona, and P. Burlamac- 
chi, IEEE J. Quantum Electron. QE-2 , 341 (1966). 

*Yu. N. Zaitsev and D. P. Stepanov, Zh. Eksperim. i 
Teor. Fiz.-Pis’ma Redakt. 6, 733 (1967) [translation; 
JETP Letters 6, 209 (1967)]. 
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FURTHER MEASUREMENTS OF QUANTUM PHASE NOISE IN A He-Ne LASER 
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ABSTRACT 

Extended measurements of laser frequency fluctuations, using 
improved instrumentation, have provided further verification of the 
Schawlow-Townes relation for the spontaneous-emission-limited line- 
width of a laser oscillator,, 
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FURTHER MEASUREMENTS OF QUANTUM PHASE NOISE IN A He-Ne Lj^ER* 

by 

R. Arrathoon and A. E-. Siegman 

Department of Electrical Engineering 
Stanford University 
, Stanford, California 


We recently- reported the observation of spontaneous-emission- 

induced phase fluctuations, or "quantum phase noise", in a He-Ne 

6328 laser operating at power levels in the low microwatt range. ^ 

Using slightly modified instrumentation, we have now extended these 

quantum phase fluctuation measurements over a wider range of laser 

powers, and have used amplitude fluctuation measurements to determine 

the excess noise factor appropriate to our* lasers. Our new results 

in general validate our earlier measurements, and give further support 

to the Schawlow-Townes relation for the spontaneous emission contri- 

2 ‘ 

bution to the linewidth of -a laser oscillator. 

The Schawlow-Townes formula predicts that the quantum phase noise 
contribution to the oscillation linewidth of a laser is 


jthf(Af ) 

<1 


N„ 


^^Ng - ^1/ 


= a 


rthf(Af )' 
' cav 


( 1 ) 


where is the full lorentzian linewidth contribution due to this 

noise source; Ng and N^ are the upper and lower laser level popula- 
tions; P is the laser power level; f is the oscillation frequency; 

-and Af is the "cold" laser cavity bandwidth. The excess noise 
cav 
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factor CC = Ng/I^Ng- (gg/g^)]?^] may be determined by measuring the ampli- 

3 

tude fluctuations in the same laser, as shown by Freed and Haus. The 
audio- frequency amplitude fluctuation spectrum of the laser, measured 
using a photomultiplier with an appropriate load resistance , may 

be written as- the sum of a shot-noise term and an excess or quantum 
amplitude noise term 


S(f) = S + S^{f). 


4 




cav 


2MTeI + 4MeI tiO , , „ 

» * ' Af y 1 + (t/Mf 


( 2 ) 


In this expression S is the spectral density of the voltage across 
2 

R (in volt /Hz); M is the photomultiplier power gain; I is the 
anode dc current; -q is the photocathode quantum efficiency; F is 
the shot noise enhancement factor; and Zif («= 10 to 30 kHz) is the 
bandwidth of the amplitude fluctuation spectrum. For lasers operating 
at low oscillation power levels, the excess amplitude noise term S^(f) 
is substantially larger than the shot noise term for frequencies less 
than Af . ‘ 

The excess amplitude noise spectrum S^(f) was measured for our 

lasers, using an audio wave analyzer, at frequencies up to approximately 

1 MHz. -The results were in general agreement with those of Freed and 
3 

Haus, except for a small additional increase in the amplitude fluctua- 
tions at very low frequencies (< 3 kHz), which we tentatively attribute 

f _ 

2 

to plasma disturbances. The mean-square voltage readings v^ = S^B , 
where B is the analyzer bandwidth, were multiplied by the appropriate 
correction factor of 1.13 to account for gaussian noise in a sinusoidally 
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calibrated .average-reading voltmeter. The value of the excess noise 
spectrum for frequencies- less than. Af , together with the measured 
■bandwidth Af of the excess noise fluctuations, in essence measures 
the excess noise factor Ct , in fact, if the measured dc voltage V 
across the anode resistance is measixred at the same time as the 

excess noise fluctuations, then the Schawlow-Townes formula may be cast 
into the extremely simple form 




3t S Af 
e 


4 / 


(3) 


a 

in which all quantities are directly measurable. The measured value of 

this expression was used to provide the theoretical value of Af^^ 

against which to compare our frequency fluctuation measurements. Using 

this approach, it is not necessary to make separate determinations of 

the photocathode quantum efficiency or the photomultiplier gain. 

However, these quantities were also independently measured to guarantee 

that all aspects of the experiment were under proper control. 

The principal change in the frequency fluctuation instrumentation 

was the replacement of the previous commercially available 30 MEIz FM 

discriminator (RHG Electronics Inc., Model DT 3006) with a quadrature 

- 4 

detector discriminator centered at 4.5 MHz. The 4.5 MHz amplifier, 

limiter, and discriminator were assembled from Fairchild and RCA 

integrated circuits. Two stages of preamplification and limiting vrere 

constructed from RCA CA-3013 integrated circuits. A subsequent stage 

consisting of a single Fairchild pATlTE integrated circuit supplied 

further amplification and limiting. This stage, in conjunction with an 

4 

external quadrature tank, also provided the frequency discrimination. 
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The narrower bandwidth and somewhat better sensitivity of this system 
permitted quantum phase fluctuation measurements at higher laser power 
levels (lower frequency fluctuation levels) than previously possible. 

Our total experimental results, together with the Schawlow- Townes 
prediction as calibrated from the amplitude fluctuation measurements, 
are presented in Fig. 1. The triangular data points indicate results 
obtained earlier by a phase jitter technique at a specific delay time 
T = 167 nsec . The square data points indicate quantum noise’ results 
obtained with the new 4.5 MHz discriminator and i.f. system. Both of 
these results differ by approximately a factor of two from the earlier 
measurements made with the 30 MHz discriminator system. The difference 
may be due to different limiting characteristics in the two discrimi- 
nators. The Schawlow- Townes predictions, taking into account uncer- 
tainties in the evaluation, seem to be in good agreement with our 
overall results . 

-7 

The apparent sharp increase in quantum noise- below P = 2 x 10 
watts should probably, as before, be discounted. In this region the 
signal-to-noise detection ratio has decreased below 20 dB; hence, all 
of the FM measurement techniques begin to deteriorate rapidly, in 
accordance with the well-known threshold properties of FM receivers. 
Also, due to the l/p power dependence of quantum noise, we might 
expect that small fluctuations in the power stabilization loop at these 
low levels could result in an augmented noise output. 

Progress is currently being made on a high-gain He-Xe 3 -51 H laser 
system with which we hope to observe quantum phase fluctuations in 
greater detail and with improved accuracy. 
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FIG. 1 — Quantiim phase noise measurements as obtained with the 
30 MHz’ discriminator, 'the ^ .5' MHz discriminator, and 
the phase jitter instrimentation. The shaded area 
represents the theoretically expected values based on 
the evaluation of Eq. (3) • 
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BUREEHT PUSHING OF 0?HE (^CILIATION FEEQUENQY OF. A 


6328 S He-Ne lASEE* 

R. Arrathoon and A. E. Siegman 

Stanford University 
Stanford, California 9^305 


ABSTRACT 

Increasing the discharge current through a small dc-excited 6328 S. 

laser chsinges the index of refraction of the discharge and causes a 

blue shift of the laser oscillation frequency of ~ 1 MHz/ma . The 

shift is independent of laser cavity tuning and decreases with increasing 

dc discharge current. Approximate calculations plus other experiments 

indicate that this shift represents the dispersive effects of the strong 

633 ^ 6402 S upward transitions from the neon- Is^ metastable 

level. The residual spectral width and FM noise in stable 6328 S. lasers, 

( 

‘especially as compared to 1 . 15 n lasers, may be due to plasma disturbances 
acting through this dispersive mechanism, rather than to microphonic 
disturbances as commonly suggested. 
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CUBEEKT FDBHIKG OF THE CSCILLA.TIOH FEEQUEKCY OF A 
6328 S. He-He LASER 

As shown in Fig. l(a), we heterodyne together the outputs from 

two gTnftll stable dc-excited. He-Ne lasers (Spectra Physics Model II9) 

12 ' ■ - 

at a nominal beat frequency of 3O MHz. ^ A small sinusoidal audio- 
frequency modulation AI added to the dc discharge current Iq of 
one laser produces a corresponding small ac frequency shift . Af in 
the oscillation frequency of that laser. This shift is' measured by 
passing the' beat note into an rf discriminator centered at 30 and 
observing the video output from the discriminator. 

Figure 1(b) plots the resulting frequency- current pushing- figure 
Af/AI" for several different values of dc current as a function 

of the tuning of the laser center frequency within the atomic linewidth, 
as measured both at 6328 S-and (by changing mirrors on the same lasers) 
at 1 . 1 %,. The 6328 5 current pushing in particular is sizeable;- decreases 
with increasing dc current ; and is constant to within measurement 

accuracy independent of laser cavity tuning within the atomic linewidth. 
The last observation indicates that the pushing represents a change in 
the background index of refraction of the laser discharge^ rather than 
any line shifts or pulling effects associated with the 6328 S laser 
atomic transition itself. 

We 'believe that this shift can be attributed to the dispersive 
(i. e., off-resonant reactive) effects of the 633^ S and 6402 R. absorptive 
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Because the metastable density saturates with increasing current, the 
slope Af/Al should decrease with increasing Iq as observed. 

As a partial test of this explanation, we focused ultraviolet 
radiation in the 3000-^000 S band from, a Hg arc lamp through one end 
mirror into the bore of one of the 6328 S lasers during operation. This 
radiation excites atoms from the neon Is up to the neon 3p levels,^ 
As expected we observed a small but measurable downward frequency shift 
in the laser frequency when the ultraviolet radiation was pulsed on, 
independent of both timing and of dc discharge current. We interpret 
the latter result as meaning that the incident ultraviolet radiation 


is almost totally absorbed on the ls-3p transitions independently 
of the Is metastable density, so that the frequency shift depends 
only on the number of incident ultraviolet photons and not on the '.Is 
density or dc discharge current , 

The observed current pushing figures for l,l 5 |i operation in Fig. 
1(b) are seen to be an order of magnitude smaller, and of opposite sign 
from the 6328 S results. They also exhibit some variation across the 
atomic linewldth. We have not pursued these results in the same detail 
as the 6328 S results, but earHl'er studies of the l.l5p transition^ 
suggest that the current pulling ‘may represent a combination of 





■3 

the strong I.O831J. transition from the metastable 2 -^S^ He level; closely 

adjacent transitions within the Ne 1 . 15 ^ oscillation group; and the more 

pronovinced pulling effects of the 1. 15p transition itself. 

Javan and -co-workers, in their "wine cellar" experiments of several 

years ego,^^^ observed a beat note between free-rxinning 1 . 15 tt He-Ne 

lasers with a spectral purity on the order of a few tens of Hz. By 

contrast, the spectral purities or noise FM bandwidths obtained in beating 

ejperimehts with 6328 S lasers generally seem to be from 2 to 3 orders 

1 2 Q 

of magnitude larger. ' ' Moreover, ih our experience improvements in 
mechanical design and acoustic isolation do not seem to yield any 
Improvement in this spectral purity beyond a certain point. Under 
quiet operating conditions, the power spectral density of the instan- 
taneous frequency fluctuation in our lasers ("G^") exhibits no discrete 

* _ s t ♦ 10“13 

oscillation spikes, such as might be attributed to plasma oscillations, 

or to mechanicsil vibration modes. The frequency fluctuation spectrum 

-2 • -2 
does show a smooth f variation, essentially the same as the, f 

noise spectrum' we observe in the plasma discharge current itself. (How- 
ever, the observed discharge ciirrent noise spectrum, when multiplied 
by the measured frequency- current pushing figure, is not sufficient to . 
account for the full observed spectral width of the laser, although both 
have the same spectral dependence). 

V 7 e suggest, ^therefore, that the residual spectral width of present- 
day high quality dc-excited (and perhaps also rf-excited) 6328 S lasers 
may be due to the dispersive effects of random fluctuations in the neon 
Is metastable density, caused by plasma disturbances and fluctuations.^^ 
These plasma fluctuations (striations, etc. ) are not presently understood 



_2 

in detail^ but In our_la$ers have an f . spectxum-and are apparently 

larger inside the discharge than is indicated by external discharge . 
current fluctuations. Substantially greater spectral purity may be 
obtainable, (with care.') in lasers^ because the discharge is much 

less optically polarizable at l» 15 ii; and hence the laser frequency is 
much less perturbed by plasma disttirbances. 
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FIG. 2 — Helim-neon energy level diagram, showing the transitions 
believed responsible for the observed frequency shifts. 
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.POSITIVE COLUMN POPULATION CALCULATIONS FOR THE EVALUATION 
. OF DISPERSIVE EFFECTS IN He~Ne LASERS 

Ro Arrathoon 

Microwave Laboratory 
Stanford University 
Stanford^ California 9^305 


ABSTRACT 

Sizeable shifts in the oscillation freq.a^ncy of a He-Ne laser with 

variations in the dc discharge current have be'^^-pb served. The shifting 

effects are believed to be due principally to changes in the populations 
3 

of the helium 2 S level and the neon Is^ level. A simplified model is 
proposed for the approximate evaluation of these populations. A slightly 
modified version of the models under regenerative conditions, is then 
used to evaluate the cizrrent shifting effects and the results are com- 
pared with experimental observations. The approach used is of sufficient 
generality to be applicable to most He-Ne lasers operat:ing under optimum 
gain conditions. The results suggest that laser operation at 1.15 P? 
rather than 6328 X, may considerably reduce the effects of current varia- 
tions on the oscillation frequency. For constant pressixre-diameter 
product conditions, theory predicts that the current- shifting effects 
will vary approximately inversely with tube diameter. 
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POSITIVE .COLUMN POPULATION CALCULATIONS FOR THE EVALUATION 
OF DISPERSIVE EFFECTS IN He-Ne LASERS 

I. INTRODUCTION 

The refraction index of a gas can he appreciably altered by the 
presence of significant excited state populations. In a gas laser 
this effect may be substantial for oscillation wavelengths in the 
vicinity of an atomic or molecular resonance. Changes in the refractive 
index of the plasma at the laser wavelength may then be effected by 
fluctuations in the populations of nearby resonance levels. The net 
■result would be to modify the effective interferometer spacing and 
produce shifts in the laser oscillation frequency. 

One method of affecting the level populations is to vary the dis- 
charge current. Earlier experiments^ at 6328 under single mode con- 
ditions, showed that low frequency perturbations in the dc discharge 
current of a small bore He-Ne laser produced blue shifts of the order 
of 1 MHz/mA. These shifts were essentially constant as a function of 
oscillation position within the Doppler broadened atomic line. Any 
pushing or p;illing effects associated with the laser transition itself 
may be expected to be markedly dependent on position within the atomic 
line. These effects should approach zero at the line center and change 
sign on either side. Since such a variation was not observed, the 
shifting was thought to be due to background dispersive effects rather 
than to any'' changes in the inhomogene ously broadened atomic line itself; 
however; an unlikely combination of asymmetric line broadening and shift- 
ing could also have produced the observed qualitative characteristics . 
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As a test of the former hypothesis, population calculations based on 
2 3 

experimental data ’ and adjusted for the conditions of our particular 

discharge were performed. These calculations indicated that background 

dispersive changes at 6328 of the proper magnitude and sign to account 

for the observed effects, could indeed be produced by changes in the 

population of the Ne Is^ level acting 'through the dispersive effect of 

the 633^ ^ and 6402 % absorptive transitions . 

For laser operation at I.I5 p, excitation induced red shifts -were 

observed that were considerably reduced in magnitude from the 6328 % 

effects; moreover, these shifts depended significantly on oscillation 

position within the Doppler broadened line. The shifting effects did 

not become zero at the line center, indicating the probable presence of 

residual background dispersive effects. Very approximate calculations, 

4 

again based on adjusted experimental data, indicated that the background 

effects were probably incurred by changes in the population of the He 

3 

2 S level acting through the I.08 p absorptive transition. 

The proposed mechanisms for the frequency shifting effects appeared 
te be consistent with the experimental observations^ and with the very 
approximate population calculations that were made. Such calculations, 
however, were not very satisfactory as large adjustments in the avail- 
able data were necessary in order to fit our particular system. Further, 
the adjustments were made on the basis of the ambipolar diffusion theory 
of the positive colimnn, which has only limited validity for most He-Ne 
lasers, -A more general and more accurate set of population calculations 
applicable to a large variety of He-He systems was clearly desirable, 

A precise calculation of the populations in a He-He system has not 
as yet been obtained and would certainly involve the simultaneous solution 



of a considerable number' of rate equations. The problem -woiold be further 
complicated by the fact that an explicit evaluation of all the terms in 
each of these equations -would probably be very approximate. An alternate 
approach is to choose a simplified model that has the proper qualitative 
characteristics and which may be expected to be approximately valid over 
most of the region of operation. .This latter approach has been chosen 
here in that it may be used to yield an approximate idea of the popula- 
tions involved, for most He-Ne lasers, with a minimum of calculational 
effort . 

The following section of this paper formulates the pertinent rate 
equations within the framework of a model that exhibits some of the 
characteristics that have been observed for He-ETe laser systems. These 
equations are then used in Section III to determine the level populations 
under both regenerative and non- regenerative conditions. The current- 
pushing effects are evaluated in Section IV and the results are compared 
with experimental data. Section V contains a- discussion of some of the 
implications of this particular -formulation. 

II. FOEMULATIOKT OF THE He- He RATE EQUATIONS 

The first qualitative description of the populations of some of 

’ 4 5 

the He-Ne levels was made by White and Gordon, ^ who found their model 
provided good agreement with the gain characteristics of the 3«39 P- 
transition. -Any qualitative evaluation of these populations requires a 
much more specific -model. In this section the plasma processes that 
must be considered in the rate equations will be explicitly defined and 
the simultaneous solution of a limited number of these equations will 
be used to describe the behavior of the relevant levels. The rate 
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equations which, shall he developed here will be derived under non- 
regenerative small signal conditions such that d;imulated emission terms 
may be neglected. 

In order to justify the simplified model that will shortly be pro- 
posedj a brief description of some of the salient characteristics of the 
level populations, the electron number densities and the electron 
temperatures in a typical He-Ne discharge will be presented. Measure- 
ments by Labuda and Gordon^ indicated that the electron temperature 
could be determined solely by the pressure-diameter (pd) product of the 
tube and that this temperature was current independent over the region 
of interest. Their measurements may be used to estimate the average 
electron temperature, for lasers in the optimum gain region of pd =3“^ 
torr-mm,^^"^ as approximately 80, 000 — 100,000°K. By' contrast, a pure 
neon discharge with typical pressures and diameters will have a con- 
siderably lower average electron temperature that might be in the 
vicinity of 20,000 — il-0,000°K.^^^ As discussed in Section III, the 
■relatively high electron temperature in a He-Ne discharge yields some 
inherent calculational advantages for the determination of level popu- 
lations. Labuda and Gordon^ also found that the electron number density 
varied linearly with current. The linear dependence of the electron 
number density on the discharge current, I , provides some justifi- 
cation for the experimentally determined relation proposed by Gordon 
k 3 ’1 

and White ■’ that describes the He 2 S metastable level population. 

This relation takes the form -Al/ (BI + C) , where A, B and C are 
constants. At higher currents the current dependent de-excitation term, 
BI , will exceed the current independent term, C , and saturation 
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will ensue. The excitation and de-excitation of the He 2^S metastahle 
level were assumed to he primarily due to collisions with electrons 
that involve only single step processes proportional to the first power 
of the electron number density. The explicit electron collisional de- 
excitation mechanism^ however, was not defined. 

There are three relevant two-body collision processes involving 
electrons and excited or unexcited atoms: excitation from the ground 

state to an excited state, excitation from an excited state to a still 
higher excited state (stepwise ■ excitation) and ionization from an 
excited state (stepwise ionization). Each of these processes has a 
converse counterpart that may be related to the forward process. The 
current dependent rate of depopulation of the He metastable levels will 
then be determined by collisional de-excitation to the ground state, 
stepwise ionization out of the excited state and stepwise excitation out 
of the excited state. Collisional de-excitation rates, which may be 
determined from the collisional excitation rates of the next section, 
are negligible in comparison with current independent loss mechanisms 
such as diffusion. Stepwise ionization rates, though somewhat larger 
than collisional de-excitation rates, also appear to be too small to 
account for the observed saturation characteristics. Stepwise excita- 
tion rates out of one excited state into another nearby excited state 
can, however, be considerably larger than stepwise ionization processes. 
If there exists a significant loss mechanism out of the nearby state 
such that most of the atoms lost to the nearby state do not return to 
the original excited state, stepwise excitation terms must also be 
included in the rate equations. 
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An energy level diagram of the helium atom is presented in Fig" 1. 

The lower four singlet and triplet levels are substantially removed from 

the next group of higher levels, and direct electron excitation rates from 

8 

the ground state to these higher levels will be significantly smaller; 
consequently, we will assume that only these lower four energy" levels 
need be considered. Helium atoms are excited to the 2^3 level by 
electron collisions and are removed from this level by stepwise excita- 
tion. There is also destruction of the 2^3 atoms by diffusion and by. 
resonant transfer collisions with neon ground state atoms that are then 
excited to the 3s^ level. Associated with this latter proves s is its 
converse, production of 2^3 atoms by resonant transfer collisions of 
ETe 3Sg atoms with helium ground state atoms. 

For stepwise excitation to be an important process in the net depopu- 
1 

lation of the He 2 3 level it is necessary to assume that some signi- 
ficant loss mechanism exists from the 2^P level so that atoms thus excited 
do not return to the 2 3 level via 2.06 p radiative relaxation. One 
possibility for such a mechanism appears to be direct relaxation of bhe 

2^P level to the ground state. The 2 ,06 inverse transition rate is in 

9 

the low microsecond range while Holstein's relations indicate that the 

1 

trapped inverse rate of the 2 P to the ground state transition is also 
in this range. Such calculations, however, are only generally indicative, 
as Holstein's results are not quite applicable to the laser situation. 

If there is trapping of the 2 .06 p transition, radiative relaxation to 
the ground state from the 2^P level will become even more important. 
Another loss mechanism for the He 2^P level appears to be resonant 
transfer collisions involving higher neon levels, several of which are 
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nearly coincident with the He 2^P level. Collisions of this kind are 
believed to be responsible for promoting laser action between the neon 
6p-6d transitions.^^ 

Spin exchange transitions between optically forbidden levels^ such 

13 ' 13 

as 2 S - 2 P stepwise excitation and 2 S - 2 S stepwise de- excitation^ 

must also be considered. The cross sections for these processes can be 

quite large and are particularly probable for transitions in which there 

12 

are no changes in the azimuthal quantum number. In fact enhancement 

3 1 

of the 2 S population by collisional de-excitation from the 2 S level 

13 

has been observed in current pulsing experiments and is believed to be 

one of the reasons for the increased power output of the 1=15 M- laser 

transition under pulsed operation. Experiments involving He-He lasers 

under steady state conditions, however, indicate that the populations of 

1 3 

both the 2 S and 2 S levels entail current dependent excitation and 

de-excitation terms that are only linearly dependent on current. ^ 

Any interaction between the 2S levels would necessarily involve higher 

powers of current. We must then assume that the upper level is somehow 

effectively decoupled from the lower level under steady state conditions. 

This will be the case if the optically forbidden interaction between the 

triplet and singlet S levels is not as inportant as the optically 

allowed Interaction between the S states and their respective P states. 

Alternately, a sufficient condition for the effective de-coupling of the 

2^S level wo-uld be that conditions of operation are such that the 2^S 
3 

and 2 S populations are related by a Boltzmann distribution at the 
electron tenperature (for a Maxwellian electron distribution) so that 
excitation and de-excitation processes between the two levels would cancel. 
In any case, we will base our model on the experimental observations that 
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1 3 

the 2 S level is essentially de-coupled from the 2 S level to obtain 

the following approximate rate equation 

dN(2^S) 

= 0 

dt 

= - Z^^I(2^S) + Z^^IT(3s2) - D^V^N(2^S) 

( 1 ) 

3 “1 

where CX^ is the direct electron excitation rate constant (cm «sec ), 

IT^ is the electron number density, is the ground state density, 

is the stepwise excitation rate constant, is the resonant 

transfer rate, Z' is the reverse transfer rate and D-, is the dif- 

fusion coefficient. The various processes that contribute to this rate 

equation are illustrated in the rate diagram of Fig. 2. 

A rate equation may also be written for the 3Sg level, in which 

the primary excitation mechanism is resonance transfer and the primary 

de-excitation mechanisms are reverse transfer and spontaneous emission. 

( 

15 

Interaction with the other 3s levels may be ignored so that the rate 
equation becomes 

dW(3s„) 

£_ = 0 

dt 

= Z^^W(2^S) - Z^^U(3s 2) - A^HOSg) . , (2) 


where is the inverse lifetime of the 3Sg level. In the vicinity 

of the tube center Eqs . (l) and (2) may be solved to yield 


N(2^S) = 


Aw + z, 

1 e Ir 


0 

^ir ^ \ 


2 . 40 ^ 1 ' 
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where a O'th order Bessel function, for the electron number density has 
been assumed, represents the on axis electron number density and R 

is the tube radius. Since most of the lasing action is confined to the 
central portion of the discharge tube, these results will be particularly 
useful when extended to regenerative conditions . The sum of the first 
two terms in Eq. (2) represents the net production rate at which 3sg 
atoms are formed due to the presence of 2^S helium. This rate, E^ , 
becomes 


E 


1 



N(2^S) 




3 

The 2 S level may be treated in a similar way. In this case 

3 3 

stepwise .excitation will occur between the 2 S and ,2 P levels . It 

3 

is not clear precisely how the 2 P level will be depopulated, as there 

are no allowed transitions between any of the triplets and the singlet 

ground state; moreover, radiative relaxation at 1.08 p from the 2 P 

3 

state back to the 2 S state will be stronger than in the case of the 
1 1 

2.06 p 2 P - 2 S transition. There will, however, be some atomic 

• 3 

collxsxonal de-excitation of the 2 P level, as .this level ,is separated 

from the 2^P level by only about 7 kT. In addition spin exchange 

electron collisional de-excitation between these closely spaced P states 

might be expected to be an even more significant interaction mechanism. 

As in the case of the 2^ level, there will also.be resonant transfer 

3 

collisions involving the 2 P level and higher neon levels. The re- 
maining excitation and de-excitation processes are similar to those for 

1 3 

the 2 S level so that the rate equation for the 2 S level will have 

the same form as that of Eq. (l) . 
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The resonant transfer rate into the neon 2 s • levels appears to be 
divided evenly between the 2sg and 2s^ levels. If we temporarily 
ignore any contributions by radiative cascade from higher levels and 
also any interactions among the 2s levels, the rate equation for the 
2s„ and 2s„ neon levels will be similar to Eq. (2). The simplest 
treatment is to consider the two levels as a group with one A coeffi- 
cient and one total de-excitation rate by reverse transfer. In this 

3 

case the solution obtained for the population of the He 2 S level 

will be similar to Sq. (3)" -It will be shown in the next section that 

Ag , the combined inverse lifetime of the 2sg and 2s^ levels, is 

1-2 orders of magnitude larger than , the reverse transfer rate 

out of these levels. As indicated in the rate diagram, the net effect 

is to make the resonant transfer process essentially unidirectional. 

Any reasonable contributions by radiative cascade from higher levels 

into the He 2s levels or from electron collisional excitation of neon 

ground state atoms to the 2s and 2s„ levels should then be expected 

^ 3 

3 

to have little effect on the population of the 2 S level. -A form 
•similar to Eq. (3) will thus be a valid solution for the population of 

3 

the 2 S level under the conditions Ag > . Any redistribution of 

populations involving all the- 2s levels would make the transfer process 
still more non-reciprocal so that the form of Eq. (3) with Ag > 
would’ correspond even more closely to the actual form of the solution 

3 

for the 2 S population. In obtaining relations of the form of Eq. (3) 
the currents are assumed to be sufficiently low so that only single step 
electron collisional processes need be considered. For very large cur- 
rents Eq. (3) will no longer be valid. 
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The treatment of the neon levels is con^licated by the helivmi 
induced population changes and also by the sheer number of interacting 
levels that must be considered. The lowest neon level, the metastable 


Is^ , appears to be somewhat -more tractable than many of the higher 

levels. For this level there will be production terms involving electron 

Gollisional excitation of neon ground state atoms. There will also be 

direct excitation of the higher groups of neon levels with subsequent 

radiative cascade contributions to the ls_ level. The’ direct excita- 

5 

tion rates of the higher neon levels, however, will be substantially 

12 

less than those of the Is levels so that these cascade contributions 
may be ignored. Helium atoms will also be supplying excited neon atoms 

1 3 

at a net production rate of , where only the He 2 S and 2 S 

resonant transfer contributions are considered. The excited neon atoms 

will ultimately radiatively relax to the neon Is levels, some fraction 

of which will become Is^ atoms. The net production rate of neon Is^ 

atoms due to the presence of helium will then be , where 

and are the respective fractional contributions arising from 

1 2 

the He 2 S and 2 S resonant transfer collisions. A rigorous deter- 
mination of kj^ and k^ requires a simultaneous solution of a large 
number of neon level rate equations. Alternately k^ and k^ may be 
approximately ‘determined by assuming that the discharge -current is 

sufficiently large' so that there is complete equilibrium among the 

2 

various Is levels at the electron temperature. In this limit the 
k's will be identical and approximately equal to the ratio of the 
statistical weight of the Is^ level to the total weight of the Is 
levels (5/12). Similarly complete equilibrium may be assumed for direct 
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electron excitation rates among the Is levels. In the equilibrium 


or high current limit stepwise excitation from the Is^ levels to the 
Ne 2p levels need not be considered in the Is^ level rate equation. 
Radiative relaxation from the 2p levels back to the Is levels to- 
gether with rapid equilibrium redistribution of the Is level populations 
ensured that stepwise excitation is not a significant loss mechanism. 
Current independent loss mechanisms for the Is^ level include 

diffusion and atoinic collisional de-excitation. Relative rate estimates 

17 

based on experimental data show that diffusion is clearly a more 
important loss mechanism than atomic collisional de-excitation. If no 
current dependent loss mechanisms are included in the formulation of 
the neon Is^ level rate equation the on axis population of this level 
the becomes 


W(ls^) = 


a N N + k,E., + k„E 
3 e n 11 2 2 


(5) 


where is the direct electron excitation rate constant, IT is the’ 

concentration of ground state neon atoms and is the diffusion 

coefficient for this level. The omission of current dependent de- 
excitation terms greatly facilitates the evaluation of the population 
of this level, but it is these current dependent de-excitation terms 
that are responsible for achieving thermal equilibrium among the Is 
levels in the high current limit. The exclusion of these terms, which 
would normally appear in the denominator of Eq.- (5), ^plicitly assumes 
that low current conditions are valid and prevents the predicted Is 
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level population from saturating. The first term of Eg. (5)^ the \m- 
saturated neon term^ will be particularly responsible for introducing 
large errors at higher values of current. The Is^ level population 
thus determined will still be very approximately correct in the interme- 
diate current region where both the high and low cxirrent approximations 
overlap. 


III. DETERMIMATIOU OF LEVEL POPULATIONS 


The following evaluation for our particular system will be generalized 
sufficiently so that similar calculations may be used for -pressure- 
diameter products in the optimum gain range. The approach used here 
will then be applicable in the range pd = 3-4 torr-mra for tube diameters 
of from l-6mm. All pertinent quantities will be determined from a 
knowledge of pressure, diameter, current and gas mixture. In the evaluation 
of the cross sections, the spontaneous emission rates, the electron 
temperatures and the electron number densities, experimentally available 
data will be used whenever possible. 

The values for electron temperature and electron ntunber densities 
may be obtained from the measurements of Labuda and Gordon^ but there 
will be some adjustment in this electron temperat-ure depending on the 
gas mixture. For H-Ne ratios ranging from 5 to 1 and higher the adjust- 
ment will be small and may be estimated- from the plots of Young's 

18 

calculations. There -will also be some adjustments in the electron 
number density measurements, which were made at only one pressiire- 
dlameter product (3.6 torr-mm) for tube diameters ranging from 2-6mm. 

These results may be extended by use of the proportionality 


•pd I 



e 


( 6 ) 
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■where is the 'electron temperature and I is the tube current. This 

relation is based on ambipolar diffusion theory and assximes that^ in the 
region of interest j the electron, drift velocity is linearly proportional 
to E/p , the ratio of the electric field to the pressure. For small 
tube diameters Labuda and Gordon-’ s measurements yield a l/d dependence 
on electron number density j however, for larger diameters the dependence 
is somewhat slower, indicating the shortcomings of ambipolar diffusion 
theory in this situation. v In the optimum gain range the tube radius is 
several times the electron mean free path, as opposed to the conditions 
of ambipolar diffusion -which require a tube radius that is very much 
larger than the electron mean free path. Equation 6, however, appears 
to have reasonable validity pro-vided it is not used over very extended 

-3 20 

ranges. For our particular tube with .a He : He ratio of 9=1 at a 
total of 3»T torr and a tube diameter of 1 mm, an electron temperature 
of 91,000° K and an on axis electron number density of 5.8-10^^ (cm"^mA.'*^) 
■was obtained. This tube is normally operated at a current of about 
5 mA. 

The electron collisional excitation and de-excitation rates may 
be obtained by appropriately integrating the electron energy distribution 
over the energy dependent cross sections of the process in question. 

The integrals may be expressed as 


OH 

e 



e m 



ct(e)E^/^ 


f(E) dE 


( 7 ) 


where (j(E) is the cross section, m is the electron mass and f(E) 
is the normalized electron energy distribution. The choice of energy 
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distribution is made particularly convenient by the relatively high’ 

average electron energy in the He-Ne discharge. For our discharge 

conditions the central portion of the distribution is indeed Maxwellian, 

but significant departures from a Maxwellian distribution are present 

21 . 

for electron energies above 35 ev. The cross sections under con- 
sideration peak -well below. 35 ev. and the electron temperature is 
sufficiently high so that most of the contribution to the rate integral 
is in the region before large departures from a Maxwellian occur. We 
will thus use the Max^^ellian as a reasonable approximation to the actual 

energy distribution in the evaluation of the rate integrals, 

s 

The direct electron excitation cross sections may be approximated 
by an expression of the form xe , which is known to have the desired 
functional characteristics. In order to obtain a better fit with the 
available data, we will use a somewhat modified version of this cross 
section. The excitation cross section will then have the form 



where is the maximum value of the excitation cross section and 

is the threshold excitation energy. The -above expression has a 
maximum at E = (C^ + C )E, and plso has the convenient characteristics 
that it may be integrated in eI\.osed form with a Maxwellian distribution. 
The cross section may then be substituted in Eq. (7) so that the resulting 
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excitation rate "becomes 



N 


e 


- Cg) 



+ 03^(2 - Cg) 



+ 20^3 


( 9 ) 

where E, , kT are in ev- The cross section for the electron 
1 ’ e 

excitation of the He 2^S level may he, determined from the partly 
theoretical^ partly experimental, compilation of Corrigan and von 

Q 

Engel. For this level the values of the parameters in the simplified 

form of the .excitation cross sections are E = 19*8 ev. , ~ 

Cg ~ .85 and C^ ~ .2 , with = 4. 5itaQ*10~^ (a^ is the Bohr radius). 

For neon the location and magnitude of the maximum value of the total 

Is level excitation cross section may he obtained from the appropriately 
8 12 

scaled data of Maier-Leihnitz. If we assume that the higher energy 

falloff is sim i lar to that of the helium 2 S level, we obtain E^ = I 6 . 6 ev. , 

2 -2 

~ .7 , *^2 — ^3 ^ ^ with Qq ~ SixsiQ -10 , The actual 

cross section of the Ne Is^ level is taken to he 5/12 of this value, 

1 

consistent with high current equilibrium considerations. The He 2 S 

level has different excitation characteristics with a much broader 

0 

maximum occuring at somewhat higher electron energies. From the 
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experimental data of Dorrestein^^ = 20.6 ev ^ ~ ~ 1 aii'i 

=i -5 ; with Oq ~ l-Sria^^ .10"^ . 

The cross sections for stepwise excitation of optically allowed 
S ■ to P states in helium do not appear to be experimentally available. 
Quantum mechanical distorted wave calculations have been derived for - such 

24 

transitions that may be approximated by a simplified relation similar 
to the classical Thompson cross section for ionization. 



where E is the energy separation between excited states, f is 
n i-K 

the oscillator strength and is the ionization potential of the 

hydrogen atom. The oscillator strengths for the 2"^ - 2^P (E^ =--.602 ev. ) 
■and the 2^3 - 2^P (E^ = I.l4 ev ) transitions are .376 and .539 respec- 
tively.^^ This cross section may be substituted in Eq,. ( 7 ) and the 
resulting expression for the stepwise excitation rate becomes 



( 11 ) 


where - E^(- is the exponential integral and E^, kT^, are 

in ev. Ihe Maxwellian velocity distributiqn assumption is more accurate 
in this case because the relatively low. value of ’E^ will make the 
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central portion of the energy distribution more important. 

The diffusion coefficients for 2 S and 2'^S helium in un- 


excited helium may be obrained from the measTxrements of Phelps 


26 


Phelps also measured the diffusion coefficients of Is^ neon in vinexcited 

5 

27 

helium and ls_ neon in unexcited neon. The values thus obtained 

5 

were -52 (cm sec .^3, .82 and .20 respectively, at 300 K and 760 
torr. These coefficients are inversely proportional to pressure, but 
their temperature dependence will be more complicated. The diffusion of 
excited atoms in a gas is somewhat smaller than that of unexcited atoms- 
in the gas due to the longer range interaction potential. Based on the 
Leonard-Jones potential considrable success has been achieved in the 
accurate calculation of diffusion coefficients of unexcited atoms over 

28 

a wide range of temperatures. The temperature dependence of the 
diffusion coefficient for excited atoms may be estimated from dimensional 
considerations that determine the temperature variation once the potential 
Interaction for a given temperature is specified. ' In the region of 
interest the potential interaction may be evaluated from Hir schf elder ’s 

28 

tabulated integrals and a slightly longer range of interaction potential 
may then be assumed for the diffusion of excited states. At constant 
density these considerations yield a temperature dependence for the 
diffusion coefficient that varies approximately as T^'^^ for both exlcted 
hellimi and excited neon in the temperature range 300-500°K. The in- 


verse diffusion coefficient for the helium and neon mixtures considered 
may be obtained by linearly weighting the individual inverse diffusion 
coefficients by their respective fractional partial pressures. Bata 
does not appear to be available for the diffusion of excited ’helium 
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in neoii but the actual value for this diffusion coefficient will be 
! 

* 2 -1 
slightly smaller than the value of 1.06 (cm sec ) for pure helium in 

27 

neon. For the particular high ratio He-We mixture considered here 

this contribution may be ignored. 

The resonant transfer destruction cross sections for the 2^S and 

2 S heliimn states involving collisions with neutral neon atoms have 

29 

been measured experimentally. The total destructive rates by resonant 
transfer, at 400°K, for the He 2^S and He- 2^S states are 


Zir = S.l^.lO^p,^ 


( 12 ) 


.. 6 * 10 ^ 


= l.b'lO" p^^ , 

2r -^Ne ’ 


( 13 ) 


where pj^^ refers to the original filling pressure at '^00 K. The 
reverse rates may be related to the forward rates and are given by 


Z' = 9.4-10 p„ 

Ir -^He 


(A) 


> 

where there is some uncertainty in the reverse rate for the second 


( 15 ) 


process due to the Interactions between the 2s and 2s levels. The 

^ J 

\ 

lifetime of the Ue 3^2 level may be determined in part by experimental 

30 31 

results and in part by theoretical calculations.-^ These results 

+7 

yield = 1. 3*10 as the inverse lifetime of this level. As measured 

32 

by Bennett-^ the inverse lifetime of the Ne 2s_ and 2s levels are 

i 

7 7 

1.0*10 aiid .63*10 respectively. These lifetimes are 1-2 orders of 
magnitude greater than the reverse transfer rate, as determined from 
equation 15* In this situation reverse resonant transfer can clearly 
be neglected with respect to the A coefficients. 
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The parameters thus determined may be substituted- in Eq. ( 3 ) to 
obtain the nonregeneratlve He 2 S and 2-^S populations. The He Is^ 
population may similarly be obtained from Eq. ( 5 ). For our .particular 
system the populations of these levels are presented in Table 1. 

Under regenerative conditions stimulated emission terms will . 
become the dominant radiative loss mechanism. The net effect will be 
to replace the term in Eq. ( 3 ) by its much larger stimulated 

emission counterpart. In the limit the term in Eqs. ( 3 ) 

and (h) will approach unity and the resonant transfer process will 
become unidirectional. The result of regenerative action will be to 
increase the net transfer rate into the He 3Sg level and decrease the 
population of both the 3^^ level and the He 2 S level. The He 2^S 
level will be essentially unaffected by I.I 5 p regenerative operation 
as the term A^/Cz^^ + A^) is approximately iinity in any case. The 
regenerative pop\ilations are also presented in Table 1. 

IV. EVALUATION OF FIffiQUENCY SHIFTIHG 


The current dependent refractive index changes at a given wave- 


length may be determined from the relation 


33 


Sn 

SI 


( 16 ) 


where \ is the wavelength^ f ^^ is the oscillator strength of the 
transition and H^ is the population of the lower level. In ,this 
expression the population of the upper level is assumed to be negligible 
in comparison with that of the lower level. ' 
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For the 1-15 transition the current dependent refractive index 
changes due to the 1.08 n absorptive transition may be evaluated from 
Eq. (l6) and substituted into the usual Fabry- Perot relation to obtain 
the following shift in the interferometer frequency 

df Sn(2^S ) 

— = - 6o6-10"' (17) 

dl Si 

In this expression we have assumed that the plasma -occupies the entire 
length of the cavity; this is nearly the case for our particular geometry. 
The current dependent population changes may be obtained from Table 1 
and the results may be substituted in Eq. (17) to obtain the values 
■plotted in Fig. 3* The experimental results are also indicated in Fig. 3 
and the two curves appear to be surprisingly good agreement^ considering 
the approximate nature of these calculations. 

In a similar way the refractive index changes at 6328 S due to the 
dispersive effects of the 633^ S and 6402 S absorptive transitions, may 
be evaluated to obtain 


Sf g N(ls ) 

— = 5,1-10 “ 2_ (18) 

6f 5l 


where minor contributions fi^^om the other’ ls-2p transitions have been 

V 

ignored; these contributions, in any case, cancel if equilibrium 
p opulat ions \ are ass\wied. 

The results of Table 1 (regenerative) together with Eq. (I8) may 
then be used to obtain the theoretical curves of Fig, 4. The curves 
indicate that the excitation dependence of the Is^ level population 
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FIG. k — Excitation induced frequency shifts at 6328 S .due to the Ke 
633k S and 6402 R transitions. The unsaturated neon value 
refers to the shifting effect in the absence of helium, 
resonant transfer contributions to the Ne Is^ level population. 
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and of the subsequent shifting effects is primarily determined by th^ 
helium in the He-Ue mixture. The experimental agreement ^ith theory 
is somewhat poorer than that obtained at 1,15 [i I however, ’ qualitative 
considerations indicate.'that if other cascade contributions to the ETe 
Is^ level population are considered and if Is^ level sat-uration terms 
are included, the agreement will prove considerably. 

Estimates of accuracy of the calculations for the He 2^ population 
evaluation can be made by utilizing these techniques for comparison with 
other experimental data. ^ On this basis we conservatively estimate 
that accuracies within a factor of two can be obtained. Our own results 
indicate that the accuracy of the He 2^S population calculations is 
probably at least this good. The He Is^ level population is not as 
accurate, partly because of the formulation of Eq. (5) and partly 
because any determination of the opulation of this level is based on 
the previous determination of helium populations. 

The experimental results, in conjunction with the population 
calculations, indicate that off resonance dispersive effects can he 
very important in He-Efe lasers. These effects are particularly significant 
at 6328 S . If the constant pd sealing consideration^ 

H I 

— « — = constant (19) 

p. d 

is substituted in Eq. (3) together with the fact that both the diffusion 
rate and the resonant transfer rates are proportional to pressure at 
constant pd, it can he seen that the population of the He 2^S level 
and the pushing effect at l,15p, is inversely proportional to the tube 
radius. Similarly E^ , and in Eq, (5) will be proportional 

to l/d , so that the pushing effect at 6328 S will also vary inversely 
with tube diameter. 
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V. CONCLUSIONS 


If the lower state population is ignored the gain for Noth the 

3.39 [0. and 6328 S transitions may be estimated from the calculated 33^ 

populations. Gains thus obtained are only of order' of magnitude accuracy 

and are too approximate for practical application. It appears that a 

7 

determination of gain based on simple empirical relationships is more 

convenient and more accurate at this time. The proposed model tends to 

overestimate the 33^ population because the assumed Maxwellian velocity 

distribution does not decrease sharply enough at higher energies. The 
♦ 

use* of a more realistic velocity distribution^ particiilarly for the 

He 2 S level which has a slowly peaking excitation cross section, should 

considerable increase the usefulness of this model. 

Order of magnitude calculations based on the Is^ level population 
indicate that trapping of radiation from the 2pj^ level to the Is levels 
is negligible. Since both the I.15 b transition and the 6328 S transition 
terminate on the 2pj^ level this indicates that radiation trapping is 
not responsible for increasing the lower state population of the lasing 
transition. 

The e^ccitation induced shifting effects were most pronounced at 
6328 £ . Residual plasma striations or perturbations in the discharge 
acting through off resonance dispersive mechanisms, may be responsible 
for limiting the spectral purity of presently available 6328 S He-Ne 
lasers. For constant pressure-diameter product conditions such plasma 
effects will be substantially reduced by the use of larger bore tubes 
operated at higher currents. 
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■Abstract — An optical heterodyne receiver is, in effect, both a 
receiver and an antenna. As an antenna it has an effective aperture or 
capture cross section ArrCn) for plane wave signals arriving from any 
direction £1. The wavefront alignment between signal and local- 
oscillator (LO) beams required for effective optical heterodyning may 
be summarized in the “antenna theorem” //An(£i)d£i= 
where the moments of the quantum efficiency t/ are evaluated over the 
photosensitive surface. Thus, an optical heterodyne having effective 
aperture An for signals arriving within a single main antenna lobe or 
Held of view of solid angle Uk is limited by the constraint 
Optical elements placed in the signal and/or LO beam paths can vary 
the trade-off between A/t and fla but cannot change their product. 

It is also noted that an optical heterodyne is an insensitive detector 
for thermal radiation, since a thermal source filling the receiver’s 
field of view must have a temperature [in (1-f^)]”' hf/k to be 
detected with S/N=l. Optical heterodyning can be useful in practical 
situations, however, for detecting Doppler shifts in coherent light 
scattered by liquids, gases, or small particles. Another antenna theorem 
applicable to this problem says that in a scattering experiment the 
received power will be ^ iVirX/djr times the transmitted power, where 
IV is the density of scatterers and a is the total scattering cross section 
of a single scatterer. The equality sign is obtained only when a single 
aperture serves as both transmitting and receiving aperture, or when 
two separate apertures are optimally focused at short range onto a 
common volume. 

Introduction 

T he fundamental virtues and weaknesses of 
optical heterodyning as a coherent optical detec- 
tion method are by now fairly well understood [1 ], 
[ 2 ]. The signal and noise properties of optical heterodyne 
and homodyne receivers have been debated [3], [4] and 
reasonably well verified experimentally [5]. The necessity 
for a stable single-frequency local oscillator has been 
pointed out [6]. The stringent alignment tolerances neces- 
sary to keep signal and local-oscillator wavefronts in 
phase over the photosensitive surface have also been 
pointed out [2], [7], [8], [9], and prescriptions for ap- 
parently relaxing these alignment tolerances have recently 
been offered [lO], [ll]. 

In discussing the directional characteristics of optical 
heterodyne detection, there is a simple and useful point 
of view based on antenna theory [2], [12] which does not 
seem to be as widely appreciated as it might be. There- 
fore, one main purpose of this paper is to point out that 
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any optical heterodyne receiver is, in essence, both a 
receiver and an antenna. As an antenna, it can be char- 
acterized by the effective aperture or capture cross section 
that it presents to a signal plane wave arriving from any 
specified direction. Moreover, and more important, we 
will prove that, in common with any antenna at any fre- 
quency, the effective aperture integrated over all possible 
arrival directions (all solid angle) is essentially X®, the 
signal wavelength squared. As a consequence, there is 
an inescapable inverse trade-off between the directional 
tolerance or the angular field of view of an optical hetero- 
dyne receiver, and its efTective aperture or capture area, 
with their product limited by the quantity X*. For optical 
wavelengths X^ is, of course, an uncomfortably small 
quantity. However, there does not seem to be any way 
to get around this limitation — relaxed directional align- 
ment requirements will inevitably bring a penalty of re- 
duced cross section or receiving aperture. 

In addition to this main point, we also wish to make two 
secondary points concerning optical heterodyne charac- 
teristics. One point is simply to note that optical hetero- 
dyne reception is virtually useless, or at least is very in- 
sensitive, for detecting the incoherent thermal radiation 
emitted by any purely thermal source. The optical 
heterodyne cannot detect less than about one photon per 
resolving time, in a single spatial mode. Thermal sources 
at reasonable temperatures emit much less than this, at 
least at optical frequencies. 

The other point has to do with the heterodyne recep- 
tion of light scattered by a volume distribution of scat- 
terers. This is a topic which has become of considerable 
interest recently in connection with optical heterodyne 
measurements of laser light scattered from small particles, 
gases, or liquids [13], [14], [15]. Optical heterodyning 
allows one to measure the Doppler shift and/ or 
the Doppler broadening of this scattered light, and there- 
by to measure the flow rate, flow gradient, and/ or various 
internal fluctuations such as critical opalescence in the 
scattering medium. 

The main point to be made here is that there is another 
antenna therorem (possibly a novel one [16]) that seems 
to apply in this case, which limits the received power from 
the scattering medium to approximately /if(rX/4ir times the 
transmitted power, where N is the volume density of scat- 
tering centers and or is the average scattering cross section 
of a single scatterer. This quantity appears to be a'funda- 
mental upper bound on the returned power in any ele- 
mentary scattering experiment using coherent detection. 
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Fortunately, this quantity, although small, is not so small 
as to rule out many interesting experiments, at least in 
the laboratory [20]. 

Tin. Anh,nna Theorem roR Optical Heterodyning 

As illustrated in Fig. 1, we consider a square-law photo- 
device which is illuminated by a local-oscillator (LO) wave 
of complex' scalar amplitude u,ix, y) and optical fre- 
quency 0 ) 11 , together with a signal wave of complex 
scalar amplitude y) and optical frequency wi. For 
simplicity, the reference plane z—0 is taken as some plane 
near the photodevice, where the signal and LO waves 
have already been combined. The z direction is preferably 
taken along a main receiving lobe or direction of maxi- 
mum sensitivity for the optical heterodyne. It is not at all 
necessary that the reference plane coincide with the photo- 
device surface, provided only that the variations in opti- 
cal path length from the reference plane to the photo- 
device surface are small compared to the difference or 
intei media te-frequcncy wavelength corresponding to the 
difference frequency (oji— wo). Transit time considera- 
tions in any practical photoelectric or photoconductive 
device will ensure that this condition is met. 


mi 


PHOTO 

DEVICE 


• Uotxjfle 




(LOCAL OSCILLATOR) 


7/' 


////^U|(x,ylB 

'll 


(SIGNAL) 

Fig. 1 . General schematic of an optical heterodyne receiver. 


For simplicity we also eliminate polarization effects by 
assuming that «o and u\ refer to the same polarization. In 
general, if the quantum ediciency of the photodevice is 
isotropic in polarization, then any incoming signal can be 
resolved into polarization components “parallel” (in a 
general sense) and orthogonal to the local oscillator’s 
polarization. The “parallel” component will then photo- 
mix with the local oscillator, and the orthogonal compo- 
nent will not. Still more generally, if the photoresponse 
(quantum ediciency) of the photodevice is anisotropic in 
polarization, then for any given LO polarization the signal 
will still be resolvable into components “parallel” and 
orthogonal (in a general sense) to the vector product of a 
photoresponse tensor times the LO polarization. We sup- 
pose here that uy dontains only the “parallel” polarization, 
since the orthogonal component will elicit no response. 
Note that with an anisotropic quantum ediciency two 
optical waves of orthogonal spatial polarization can beat 
together [12], [17], [18]. Also, in either isotropic or 
anisotropic cases if a signal wave is randomly polarized 
with a uniform distribution, half of the signal intensity 
will always remain undetected. 

* Any quaiUity having a superscript tilde will, in general, be 
complex. 


The scalar 'amplitudes are assumed to be normalized 
so that the incremental photocurrent, produced at the 
pholodevice’s output terminals by tiie complex scalar in- 
tensity u(x, y, t) falling on an incremental area dA on 
the reference plane, will be 

dl(l) = v(-f,u)\ !j,l)\-<{.\ 

or 

■f(0 = J'J* vU‘, !/) I «(.r, .'/, 0 I'-dA, 

where tj(^x, j) is the quantum ediciency for light striking 
the point (x, y). in tlie elementary optica! heterodyne 
case, the total scalar intensity is 

«(t, !/, I) = «a(.!‘,'.'/)c^““' + «i(.r, ,v)f-'“'>'. 

Hence the total photocurrent is 

/(i) = /.) + f\ + -f- Iv,*r 

where the dc currents /o and 4 induced by the local oscil- 
lator and the signal are given by 

/() = J'J' v(x, !j) I Ha(.r , !/) |Wzl 
= J'J' !/) j «i(.r, !j) \i(lA, 


and where the complex phasor amplitude of the difference- 
frequency photocurrent is given by 


ijJ 10 


J J '/)"i(.f. .v)«o*(r, u)dA. 


If M(i(jc, y) and Ui(x,y) have the same spatial variation (in 
magnitude and phase) over the reference plane, this re- 
duces to the optimum photomixing situation 

I = Alufy 

which gives maximum sensitivity [3], [4]. 

Suppose now that the incident signal wave is a uniform 
plane wave with wave vector ky^(k^, k„, k.), so that 

«i(.r, If, z) = 

On the reference plane this gives 
]J, 2 = 0)= 

If the signal wave vector has magnitude /ci= j Ai[ 
and if the direction of arrival of this wave is along spher- 
ical coordinates 6 and <j> in the coordinate system of Fig. 
I (where 9=0 is along the z axis and ^ is measured from 
the X axis in the x, y plane), then 

A'j. = kysiw 6 cos <f> 
k„ = ky sill 6 sill <j>. 

The difference-frequency photocurrent is then given by 

iho = n r r ’?(•'■) .'/)«»*(•<■, 
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Now we define a mean quantum efficiency ij given by* 


t] S 


f f y) 1 ‘ticix, y) I’^dxdy 

JJ \u,i{x,y)\^dxdy 


Using this definition and the previous equation, we choose 
to write the magnitude squared of the IF photocurrent in 
the form 

|Aol* = 4/o[^leihAH]. 

The product n|eil* iit this expression is the dc photo- 
current density, i.e., the dc photocurrent that would be 
produced by a plane wave of complex amplitude el falling 
on a photodevice of unit area and quantum efficiency ij. 
The whole quantity in the square brackets is the effective 
signal current that mixes with the LO dc current to pro- 
duce an IF photocurrent. This effective signal current will 
be in generi less than the actual signal dc current, which 
will be i/|el[* times the actual photocathode area. The 
quantity Ar has the dimensions of an area, and in fact 
may be interpreted as the effective receiver area or re- 
ceiver aperture of the optical heterodyne receiver in re- 
ceiving a plane wave signal by mixing with the LO wave. 

This effective receiver aperture Ar is closely analogous 
to the effective aperture of an ordinary radio wave an- 
tenna. It is a function of the vector arrival direction of 
the signal wave and, from the previous three equations, 
is given by 


Ar(JCx, kii) — 

I J*J* ij(a:, y)u*{x, y)e-‘^^’^^+’='"''>dxdy 




From the form of this expression, the effective aperture 
An^kx, kji) is proportional to the power spectral density 
in the space frequencies kx, ky of the spatial distribution 
ij(x, y)uo*(x, y), i.e., the LO complex amplitude weighted 
by the quantum efficiency. 

We now wish to evaluate the integrated value of 
over all possible directions of arrival Q. The dif- 
ferential of solid angle can be written 

da = sin 6ddd<ft. 


But we can also note that, from the relations above, we 
can write 

dkxdk,, = fci* cos 6 sin 6d8d(}> = cos 8dCl. 


* We choose this definition for i because, with this definition, the 
minimum detectable signal power with an optimum spatial variation 
of Hi(x, y) wilt be just photons per reciprocal bandwidth (based 

on S/N=l, shot-noise-limited). 


As mentioned earlier, we suppose that an optical hetero- 
dyne receiver will usually have only one (or at most a few) 
do minan t (and reasonably narrow) main receiving lobe(s) 
pointed along the 6=0 or 2 axis, so that over the range 
of an^es for which Aidpt) is appreciable we can make the 
approximations cos and dSif^^kr^dk^dky. We can 
then write 


//. 


all solid angle 


AH(o)(ia 


(lAt*) 



eo 


But, from the interpretation of Ai^k^y kji) as a spatial- 
frequency power density we can easily find that^ 






J* J* ^ \^dcdy 

fi* i uaix, y) \ Hxdy 


or 


If 


AR(ii)d£i 



where ij* is implicitly defined in the previous equation. 

If the quantum efficiency is uniform over the photo- 
device, then the ratio is strictly imity; and in any 
reasonable case this ratio is not likely to be far from unity. 
We finally conclude, therefore, that the integrated effec- 
tive aperture of an optical heterodyne receiver is given to 
a good approximation by the simple relation 


// 


A«(Q)dO « X* 


If the antenna has a single main receiving lobe that sub- 
tends a solid angular field of view of Hr steradians, with 
an effective aperture of .<4^ for sources inside this field 
of view and zero outside, then [2], [12] 

ArHr ~ X*. 

These relationships are, of course, identical with the well- 
known antenna theorem for the integrated cross section 
or integrated gain of a single-mode antenna at any fre- 
quency. A major problem is the extremely small size of 
the wavelength at optical frequencies. 


Examples 

The simplest illustration of this antenna theorem is a 
photosurface of physical area A (diameter d) uniformly il- 
luminated by a uniform plane-wave LO beam, as shown in 


* We are basically using the Wiener-Khintchine relation for 
Fourier transform pairs. 
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(c) 



Fig. 2. Several sketches illustrating various concepts 
in optical heterodyne antenna theory. 

Fig, 2(a). The effective aperture for a signal plane wave 
parallel to the LO plane wave is clearly just A/{=A. If the 
phase variation between signal and LO waves is not to 
exceed one wavelength across the aperture, the wavefront 
must not be tilted from parallelism by more than Ad!=>\/d. 
Therefore the incoming signal wave must be confined to a 
cone of solid angle given by 2h^(A9)-^\-/A. Therefore, 
we have 


^ AX « V-, 

in agreement with the theorem. 

Figure 2(b) illustrates an arrangement that might at 
first seem to overcome the antenna theorem limitation. 
The signal wave is collected and focused by a lens of area 
A;, and focal length / onto a cathode of area Ac located 
along the focal surface of the lens. The signal spot size at 
the focal plane is As. The LO beam diverges as a spherical 
wave of proper curvature so as to uniformly illuminate 
the focal-plane cathode at normal incidence. 

The matching of wavefronts in the overlap region of 
signal and LO beams can be made as ideal as desired, at 
least in principle, for a signal spot anywhere on the 


cathode surface. Therefore, the solid angular field of view 
is given by 



At the same time it would appear that the effective aper- 
ture can be made indefinitely large simply by increasing 
the lens area A The restriction in the present case comes 
about, however, because the lens area Ai, and the focal 
spot size As are interrelated. For ideal optics the diffrac- 
tion limit of the lens is 

X- As 

'a'j.'Y' 

The LO dc current la and the related shot noise are 
determined by the total LO intensity over the entire 
cathode. However, the signal light mixes only with the 
fraction As! Ac of the local-oscillator light lying within 
the focal spot. In effect, therefore, the conversion gain of 
the heterodyne process is reduced by this fraction. The 
effective aperture is, therefore, given by 

The net result is 



which again agrees with the antenna theorem. ' 

We close this portion of the discussion with three gen- 
eral remarks. 

1) First, as Fig. 2(c) illustrates, an Af-power telescope 

placed in the signal beam only, in front of an 
optical heterodyne, will magnify the apparent area of the 
cathode by M-\ but at the same time the angular sensi- 
tivity will be increased by Af-. Hence with such a telescope 

AAhi' « (il/M/0 X 

In general, optical elements of any sort placed in the sig- 
nal beam before it enters the optical heterodyne will 
change, both the effective aperture An and the angular 
field of view iht, but always subject to the basic constraint 
that 

2) Optical elements of any sort placed in the signal and 
LO beams after they are combined will in general not alter 


•* References [10] and [1 1] give an experimental demonstration of 
essentially this scheme, in conjunction with an image dissector photo- 
device that passes only the photocurrent from the focal plane spot 
As- Their final system also automatically tracks the movements of 
the spot Ah. This system does not violate the antenna theorem, 
nonetheless, because its msiwitmeous field of view is still very nar- 
row, so as to allow a relatively large aperture. The resulting critical 
alignment is made more convenient by the electronic and servo- 
controlled alignment method, but the alignment is still critical. Also 
a large waste of LO power is inherent in the scheme as used. 
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either Ar or £!« (provided that the total light intensity 
continues to strike the photosensitive element). Figure 
2(d) illustrates this for a very simple example. 

3) The integral expression given above for Adkx, in 
terms of ij(x, y)«o(x, y) has exactly the form that would 
result from using Huygen’s principle together with the 
Fraunhofer approximation to evaluate the intensity dif- 
fracted into a direction *1 from a source amplitude distri- 
bution ij{^, j)j3o(x, y). Therefore, the following is a pre- 
scription for determining the antenna pattern of an opti- 
cal heterodyne [12]. Consider the complex LO amplitude 
distribution falling on the photodevice surface (weighted by 
the quantum efficiency distribution if necessary). Reverse the 
direction of propagation of this LO distribution and allow 
the reversed wavefront to propagate back out through any 
optical elements that an incident signal wave would traverse. 
The resulting far-field or Fraunhofer diffraction pattern will 
be the antenna pattern of the optical heterodyne receiver 
[19]. 

Optical Heterodyne DETEcnoN of Scattered Light 

As noted in the Introduction, one application where 
optical heterodyning has proven very fruitful is in making 
remarkably detailed measurements of the Doppler shift 
and/or the Doppler broadening of laser light scattered 
from gases, dust particles, and other scattering elements 
[13], [14], [ 15 ]. Therefore, we will next show that there 
appears to be a very simple and general expression for 
the maximum signal return to be expected in this type of 
experiment. 

We consider an experimental situation such as that 
shown in Fig. 3, where the transmitter is a laser, and the 
receiver is an optical heterodyne with the usual field of 
view fl/i and receiving aperture Ar. The transmitter and 
receiver beams intersect in a common volume which is 
presumed to contain a density N per unit volume of scat- 
ters. Each scatterer is assumed to have a total scattering 
cross section tr; for simplicity the scattering is assumed to 
be isotropic. Only single scattering is considered. 

For simplicity we assume rectangular transmitter and 
receiver beams, as shown, so that the scattering volume 
common to the two beams is a parallelopiped with sides 
a, b, c. Then, the power density in the transmitter beam in 
the scattering region is sin fl, and the total power 

scattered by an individual scatterer is (<r/ab sin 0)Ptmns- 
Of this scattered power, a fraction Ar/4ttR^ will be inter- 
cepted by the optical heterodyne. The total number of 
scatterers in the common volume is Nabc sin 0, and the 
receiver field of view is equal to ^R—bc sin 0/R^. By com- 
bining all of these factors, we obtain 

f^rpceivcd “ (iV * ohc slu 0){ truns 

V \ab sin 9/ 

X Na\ ^ 
osmd 4jt 



Fig. 3. Geometry for the analysis of heterodyne 
detection of scattered light. 


If the transmitter and receiver are both focused on the 
scattering region to give a minimum value of b, the result 
with an / 7 I lens will be ftssX. Moreover, the factor 
(sin 6)“^ does not allow any large improvements, particu- 
larly in the focused case, because the region of beam over- 
lap will be limited by the depth of focus. Therefore, by 
setting \/b sin 0^1 we are led to the simple limiting ex- 
pression [ 20 ] 

I^rccelvcd ^ <rX 

Ftrnn. ~ 4 t 


This same result can be obtained directly by analyzing the 
situation in which transmitter and receiver are focused 
onto a common small volume. 

This conclusion is reinforced if we consider as an alter- 
native the case of a single transmitting and receiving aper- 
ture which we assume to be imbedded in the scattering 
medium. The geometry for this case is shown in Fig. 4. 
In this case the transmitted power density at a range r in 
the far field of the aperture (assuming the system to be 
focused at infinity) is Ftrnns/r*B«- The power scattered by 
one scatterer is (cr/i^ClR)Piram, and the fraction of this 
which is detected is AR/Airr^ The number of scatterers in 
a shell of thickness dr at range r is NSlnr^dr. Therefore, we 
obtain for the received power from this shell 

dPr,^i^d = (NilarMr) (—^ (— — 

\47rrV \rW 


NctAr 

4t 


dr 

“ P trans. 

J.2 


This expression holds only for ranges r in the far field 
of the aperture, beginning at about one “Rayleigh range” 
ro from the aperture. The “Rayleigh range” of an aperture 
is defined as the distance where the diffraction spread 
/"o’Qb from the aperture is approximately equal to the 
original aperture size Ar, or 


ro « 



Ar 

X 
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Fig. 4. Heterodyne detection of scattered light using the 
same aperture for transmitting and receiving. 


derivation, without altering the results if the source is 
assumed to be planar and located at a range R in the far- 
held region of the receiver. The total power AP emitted 
into- a hemisphere, within frequency range A/, from a 
small area AA on a blackbody is 


AP = 


27t hfAf 
AA. 

X2 - 1 


Now, scatterers located closer to the aperture than about 
one Rayleigh range are “out of focus.” That is, they scat- 
ter spherical waves which do not become sufficiently 
planar at the receiving aperture to be properly detected. 
Moreover, inside the Rayleigh range — that is, in the 
Fresnel region — the transmitted power density does not 
continue to rise as 1/r®, but rather varies between zero and 
the value at about one Rayleigh range. 

It seems a reasonable approximation, therefore, to 
neglect those scatterers located inside the Rayleigh range, 
and to sum up the contributions only from shells located 
at and beyond the Rayleigh range. The result is 

Ptcctnv&i ^ ^ Aji p ^ d.T AftrX 

Pit ms 4 tT «/ r=rs 47T 

This single-antenna result is the same as the limiting value 
of the two-antenna analysis just given, and is independent 
of the optical heterodyne’s aperture or field of view. The 
same result is also obtained for a single antenna focused 
at any arbitrary range in the scattering volume [20]. 

This limiting formula, derived only by these approxi- 
mate arguments, appears to be a general limitation on 
scattering experiments [l9]. The expression given is use- 
ful in evaluating the minimum scattering density^ Nor that 
can be detected with-a given transmitter power and opti- 
cal heterodyne receiver. Fortunately, there are situations 
of much practical interest [13], [14], [15] in which very 
useful scattering measurements can be made taking ad- 
vantage of the unique-capabilities of the optical hetero- 
dyne receiver for detecting Doppler shifts and Doppler 
broadening of the scattered light. 

Heterodyne Detection of Thermal Radiation 

One question that seems to arise frequently is the utility 
of the optical heterodyne receiver for detecting optical 
radiation from incoherent thermal sources. This question 
is basic, for example, in determining whether optical 
heterodyne detection has applications in astronomy or 
radiometry. We will next give a simple argument to show 
that optical heterodyne detection is a relatively insensi- 
tive detection method for an incoherent thermal source. 

Consider an optical heterodyne receiver of effective 
aperture An whose field of view Qr is filled by a black- 
body thermal source at temperature T. It simplifies the 

‘‘ In scattering from index variations in liquids, N<r corresponds 
directly to the turbidity. 


Since the power per steradian emitted normal to the sur- 
face is l/ir times this, while the solid angle subtended by 
the receiver .aperture as seen from the blackbody is 
Alt/ R?, only a fraction Ar/ttR} of the blackbody radiation 
is received by the optical heterodyne; and only ^ of that 
on the average is properly polarized so as to be detected. 
Since emission from different points on a thermal source 
is incoherent, we can simply add the power contributed by 
each small area AA within the total area R?Qr on the sur- 
face that lies within the field of view of the optical hetero- 
dyne receiver. 

The net result is that the average signal power received 
by an optical heterodyne from a thermal source filling the 
field of view is 


P ■ = 

■*- SIS 



hfdf 

gh/IkT _ 



hfdf 

gh/IlT _ I 


But, the equivalent noise input to an optical heterodyne 
receiver, due to LO shot noise, is well known to be 

P.OUO = 0/v)hfdf. 

Hence, the signal-to-noise ratio for optical heterodyne 
detection of the thermal source is 

S n 


For the thermal source to cause a signal-to-noise ratio of 
unity or better, the minimum source temperature must be 

hf 1 hf I 

T>— ij«i. 

k ln(l -j- k 1 ] 

For frequencies in the visible light range, the quantity 
hf/k is about 25 000°K. Therefore, it would require a 
source I /ij times this hot, filling the receiver field of view, 
to cause the noise power level in the heterodyne IF chan- 
nel to double. To put this in another way, an optical 
heterodyne receiver pointed directly at the sun (6000°K in 
the visible) would show virtually no additional signal. 

Detection capabilities for thermal targets can, of 
course, be considerably increased by using chopping 
and time-integrating techniques, and other radiometric 
methods. Also, the relative sensitivity or the minimum 
detectable source temperature for an optical heterodyne 
improves directly with increasing wavelength as one goes 


October 1966 / Vol. 5, No. 10 / APPLIED OPTICS 1593 



into the infrared. Nonetheless, it seems clear that the opti- 
cal heterodyne method does not offer high sensitivity for 
detecting purely thermal radiation.® 

Conclusions 

The optical heterodyne receiver is usefully viewed as 
both a receiver and an antenna. As the latter, it shares, 
with all other antennas, an unavoidable reciprocal trade- 
off between effective aperture Ar and angular field of 
view or angular alignment tolerance fix; the product 
The heterodyne receiver is relatively insensi- 
tive for detecting incoherent thermal sources. It is effec- 
tive for detecting Doppler shifts, especially in scattered 
coherent laser light. ITiere is a simple limiting expression 
which says that the scattered light received by an optical 
heterodyne cannot exceed N&X/Att times the transmitted 
light, where N is the density of single scatterers and <r is 
the total scattering cross section of a single scatterer. 
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